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We show that the idea of mapping between the Newtonian and Brownian diffu- 
sivities proposed and tested on a class of particle systems interacting via soft and 
ultra-soft potentials (IPL, Gaussian core, Hertzian, and effective star-polymer) by 
Pond et al, [Soft Matter 7, 9859 (2011)] is also applicable to the Yukawa (screened 
Coulomb) interaction. Some of the implications of this result with respect to self- 
diffusion in strongly coupled complex (dusty) plasmas are discussed. 

PACS numbers: 52.27.Lw, 66.10. C-, 66.10.cg 

I. INTRODUCTION 

Recently Pond at a/.- proposed a simple idea for mapping between long-time self-diffusion 
coefficients obtained from molecular and Brownian dynamics simulations. Their primary in- 
terest was on particles interacting via soft inverse-power-law (IPL) potentials as well as 
ultrasoft Gaussian core, Hertzian and effective star-polymer interactions. The latter three 
models have nowadays received considerable attention in the field of soft condensed mat- 
ter. A simple heuristic expression relating Newtonian and Brownian diffusivities has been 
demonstrated to describe reasonably well extensive simulation data for the class of systems 
investigated.- 

The main purpose of the present work is to check whether the proposed mapping is ade- 
quate for complex (dusty) plasmas, which have been recently recognized as new class of soft 
matter - the "plasma state of soft matter" .— ^ Complex plasmas consist of weakly ionized gas 
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(conventional plasma) and highly charged macroscopic (dust) particles.- - - Highly charged 
particles interact with each other electrically. In the first approximation this interaction 
can be modeled by the Debye-Hiickel (screened Coulomb or Yukawa) potential, although its 
actual shape can be considerably more complicated, especially at large interparticle separa- 
tions.-~- The electrical interaction energy can often be remarkably high as compared to the 
particle kinetic energy and, therefore, complex plasmas can be viewed as classical systems 
of individually visible strongly interacting particles. They can be used (complimentary to 
other soft matter systems like e.g. colloids and granular medium) to investigate a broad 
range of important fundamental processes, such as phase transitions and self-organization, 
phase separation, rheology, waves, transport, and other. 

The focus of the present study is on the diffusion-related properties of complex plasmas. 
We have to admit that diffusion in Yukawa systems has been already investigated in a number 
of papers. For instance, molecular dynamics technique has been employed in Refs. Iioljll to 
simulate self-diffusion coefficients of one-component Yukawa systems in the fluid phase in a 
wide range of thermodynamical parameters. Brownian dynamics technique has been used in 



Ref. 
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to obtain the diffusion coefficient of an overdamped Yukawa fluid near the freezing 
point. Comparison between molecular dynamics and Brownian dynamics diffusivities in 
a system of polydisperse Yukawa particles near a glass transition towards an amorphous 



solid has been performed in Ref. 
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A Langevin-type approach has been used to study 
diffusion in the fluid phase of Yukawa systems over a broad range of damping strength.— r — 
These latter studies are particularly relevant to complex plasmas since the coupling strength 
between the particles and the plasma background (i.e. damping, which is usually dominated 
by the neutral gas component) can be varied in a very broad range. As a result, complex 
plasmas can be in principle "engineered" as essentially a one component system (when 
damping is vanishingly small and the Newtonian dynamics is dominated by the interparticle 
interactions), or as an overdamped system of Brownian interacting particles (when damping 
from the background medium dominates dynamical phenomena).— This explains why the 
possibility to map between molecular dynamics (Newtonian) and overdamped (Brownian) 
diffusivities represents considerable interest in the context of complex plasmas. In the rest 
of this paper we use the available numerical results to demonstrate the applicability of the 
mapping proposed in Ref. Q to strongly coupled Yukawa systems. Moreover, we take a next 
step and suggest an expression for the diffusion coefficient, which is applicable to strongly 
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coupled Yukawa systems with an arbitrary damping strength. 

II. NORMALIZATIONS 

The Yukawa pair potential can be written as 

U{r) = e(A/r)exp(-r/A), (1) 

where e is the energy scale, A is the screening length, and r is the distance between two 
particles. For charged particles immersed in a plasma the energy scale is related to the 
particle charge Q (e = Q 2 /X) and the screening comes from the equilibrium redistribution 
of plasma electrons and ions in the vicinity of the particle (A = Ad, where Ad is the Debye 
radius). The Yukawa system in thermodynamical equilibrium can be characterized by two 
dimensionless parameters. One natural choice is the reduced temperature T* = T/e and 
density p* = pA 3 . In the field of complex (dusty) plasmas, however, it is more usual to use the 
coupling parameter T = Q 2 /(TA) and the screening parameter n = A/A, where A = p -1 / 3 
is the mean interparticle distance. Note that the fluid-solid phase transition in Yukawa 
systems can be to a good accuracy described by T ~ 106exp(K)(l + k + |k 2 ) -1 , provided k 
is not too large.— We will denote the value of the coupling parameter corresponding to this 
phase transition as T M , where the subscript "M" refers to melting (similarly, T M corresponds 
to the temperature at melting).— 

In the following we will be dealing with the diffusion coefficients of 3D Yukawa systems 
in the fluid phase only. The actual (dimensional) long-time diffusion coefficient is denoted 
by D, its values in the limiting regimes of molecular (Newtonian) dynamics and Brownian 
(overdamped) dynamics are denoted by -Dmd and -Dbd respectively, following the notation 
of Ref. Q. For Newtonian diffusion without damping we employ Rosenfeld's normalization^ 
and use the reduced diffusion coefficient in the form -D RM d = -DmdP 1 ^ 3 V M /T ', where M is 
the particle mass, and the temperature T is in energy units. For Brownian diffusion the 
natural choice of reduced diffusivity is -Dbd/ Dq, where Do is the value of -Dbd in the dilute 
limit (or, equivalently, in the absence of interparticle interactions). 
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III. MAPPING 



A heuristic expression describing mapping between Newtonian and Brownian diffusivities 



proposed in Ref. 



is 

.3/2 



D BT3 /D = 1 - (1 + c^rmd + c 2 Cd) > (2) 



where c\ = 3.3176 and C2 = 2.6645 are numerical constants. It describes relatively well 
the numerical data for IPL interactions as well as for model complex fluids with soft star- 
polymer, Gaussian-core, and Hertzian interactions.- 

In order to verify whether this mapping is applicable to Yukawa interactions, let us first 
consider the available numerical data for the Newtonian limit. Figure [T] shows the reduced 
diffusion coefficients -D RMD plotted against the ratio T/r M = T M /T, 1.6. ? clS 8b function of 
the distance from the melting line. Numerical data for Yukawa systems are taken from 
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the numerical values 



MD simulations by Ohta and Hamaguchi.— As discussed in Ref. 
of the diffusion coefficient may be subject to errors of up to about 10%. The data for the 
limiting case (k = 0) of one-component-plasma (OCP) are those from MD simulations by 
Hansen et al— Figure Q] demonstrates that the reduced diffusion coefficient -Drmd shows 
a quasi-universal behavior, except in the regime of sufficiently steep interaction (k > 4) 
and weak coupling, where systemmatic deviations (decrease in -Drmd) are evident. In the 
strongly coupled regime not too far from the melting line (say, 0.1T < T < Tm) the data 
are scattered near a single curve (see inset in Fig. [TJ), which can be reasonably fitted by the 
function 

Drmd ^ a (T/T M - if + 7, (3) 
with the coefficients a ~ 0.05, /3 ~ 0.8, and 7 ~ 0.03. The same functional form has been 



used in Ref. 
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but with different coefficients, because of the different normalization for the 
diffusivity employed there. 

With the use of Eqs. (T5])-(j3]) we can immediately calculate the dependence of the reduced 
Brownian diffusion coefficient D-qb/Dq on the ratio T /Tu (and Tu/T). The results are 
shown in Fig. [2] by the solid curve. Here we also show a simple linear dependence 

Dbd/Dq ~ 0.3(1 - T u /T) + 0.1, (4) 

proposed by Vaulina et al— for strongly coupled Yukawa systems (r > 0.5r M ) in the 
limit of strong damping (dashed line) along with the numerical data from the same work 
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FIG. 1: Reduced diffusion coefficient -Drmd in the molecular dynamics (Newtonian) regime vs. the 
relative coupling strength V/Ty\. Numerical data are shown by symbols. The data for the Yukawa 



potential are taken from Ref. 
Ref. 
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. The data for the one-component-plasma (OCP) limit are from 



. Inset shows the portion of the data corresponding to the strongly coupled regime. The 



solid curve is the fit using Eq. ([3|). 

corresponding to the highest (but finite) damping rate (circles). The agreement is convincing, 
which proves the applicability of the mapping ([2]) to Yukawa systems in the strongly coupled 



regime. 



IV. DISCUSSION 



Let us discuss an important consequence of the applicability of the mapping ([2]) to strongly 
coupled Yukawa systems. First of all note that on approaching the melting line, both -Drmd 
and -Dbd become small (see Figs. [Hand [2]). Hence, to the lowest order in .Drmd we have 

AWA) - ciAimd, ( 5 ) 

i.e. both reduced diffusivities are proportional to each other. As we have seen already, in a 
broad range of parameters corresponding to strong coupling, the value of -Drmd is governed 
by a single quantity - the ratio of r/T M (or T M /T). Let us denote J r {T/T M ) = CiDrmd- 
Then the actual diffusion coefficient in the two limiting regimes considered can be presented 
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FIG. 2: Reduced diffusion coefficient D-qd/Dq in the overdamped (Brownian) regime vs. the 



relative coupling strength T/Tm- Numerical data shown by symbols are from Ref. 



15. The solid 



curve represents the mapping by Eq. ([5]). The dotted curve corresponds to its simplified version 



). The dashed line is the linear fit suggested in Ref. 



15l for the strongly coupled limit. All three 



curves and symbols almost fall on top of each other for T •> 0.5Fm- 



as 



D = Dn 



(6) 



f -F(r/r M ) (Newtonian) 

[ .F(r/r M ), (Brownian) 
where b>f v is the macroscopic friction rate associated with particle interaction with the sur- 
rounding medium (neutral gas in complex plasmas), so that the bare Brownian diffusion 



coefficient (no interaction) is Dq = T/(Mvf I ). The "damping index" £ = u^A-^ M/T natu- 
rally measures the damping strength. Physically, £ is the ratio between a typical interparticle 
spacing and the mean ballistic free path of the particles. For ^ C 1 the particle motion 
is essentially ballistic and this corresponds to the molecular dynamics (Newtonian) limit. 
Contrary, for £ ^> 1 the particle motion is diffusive on the length scales considerably shorter 
than the interparticle distances. This corresponds to the overdamped (Brownian) limit. Ex- 
amination of Eq. ([6]) immediately suggests that it is convenient (and not very unreasonable) 
to separate the effect of interparticle interactions and damping and write the actual diffusion 
coefficient in the form 

d = D Q T(r/r u )g(t), (7) 

where the function J 7 is responsible for the suppression of self-diffusion due to (strong) 
interparticle interactions, while Q is the "damping function" describing the transition from 
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FIG. 3: The damping function £/ vs. the damping index £ = i/fyAy 7 M/T for strongly coupled 



15. The solid curve is a 



Yukawa systems. Symbols correspond to the numerical data from Ref. 
simple approximation given by Eq. ([SJ. 

the limit of Brownian dynamics [£ 3> 1 and Q (£) — > 1] to the limit of Newtonian dynamics 
[£«land 0(£)->£/ Cl ]. 



We use the numerical data from Ref. 
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to check the plausibility of the scaling (J7J). These 
data were obtained for the two values of the screening parameter (k = 2.42 and k = 4.84) 
and various damping indexes 1 ^ £ ^ 100. In figure |3] we plot the corresponding ratios 
D/lDoJ^lT/Tu)} = Q(0- Although the data points are somewhat scattered, they clearly 
tend to group near a certain curve. Thus, in the regime of strong coupling, Eq. (J7J) represents 
a reasonable approximation (we found it, however, less accurate for weaker damping). A 
simple function of the form 

0(0 = 7T7 (8) 



is exact in the respective limiting cases and provides a reasonable description of the numerical 
results in the transitional regime. More statistics and higher data accuracy are apparently 
required to construct a better analytical expression for £/(£). 



V. CONCLUDING REMARKS 



We have demonstrated that the idea of mapping between the Newtonian and Brownian 
diffusivities proposed and tested on a class of soft repulsive potentials (IPL, Gaussian core, 
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Hertzian, and effective star-polymer) in Ref. Q, also works relatively well for the Yukawa 
interaction potential. This can have important implications for the field of complex plasmas, 
where the coupling strength between the particles and surrounding background (typically 
dominated by neutral gas) can be varied in a broad range, so that both limiting regimes (as 
well as intermediate states) can be realized. 

We have also seen that in a wide parameter regime the Newtonian diffusivity -Drmd of 
Yukawa systems exhibits universal behavior with respect to the reduced coupling strength, 
provided the interaction is soft enough (re < 4). An interesting issue to investigate is whether 
this or similar scaling holds for other soft and ultrasoft interactions, for which the universal 
excess entropy scaling^ does not work.->22 

We have proposed an expression for the diffusion coefficient applicable to an arbitrary 
damping efficiency in the system. In this expression the ratio of the actual diffusion coef- 
ficient to its bare Brownian value is given by a product of the two functions, one of which 
describes suppression of the self- diffusion due to interparticle interactions and depends only 
on the relative coupling strength, while the other - the damping function - is a convenient 
measure of the damping strength and is independent of the interaction details. Such an ap- 
proach demonstrates reasonable accuracy when applied to Yukawa systems in the strongly 
coupled regime. Thus, it can be used to estimate the diffusion-related quantities in complex 
plasmas near the fluid-solid (crystallization) phase change. The important question to be 
answered in future is whether this scaling is truly universal and/or how wide is the class of 
interaction for which it holds. 

Finally, it would be interesting to look at the observed universalities from the perspective 
of the emerging principle of correspon ding states for strongly coupled particle systems, which 



has been recently put forward in Ref. 
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